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Abstract 



P^ This work is a follow-up to our previous work [2]. It extends and complements, both 

•^r theoretically and experimentally, the results presented there. Under consideration is 

the homogenization of a model of a weakly random heterogeneous material. The ma- 
terial consists of a reference periodic material randomly perturbed by another periodic 
material, so that its homogenized behavior is close to that of the reference material, 
y We consider laws for the random perturbations more general than in [2 . We prove 

the validity of an asymptotic expansion in a certain class of settings. We also extend 
I the formal approach introduced in [2 . Our perturbative approach shares common 

^ features with a defect- type theory of solid state physics. The computational efficiency 

C^ of the approach is demonstrated. 
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T! 1 Introduction 

k>< Our purpose is to follow up on our previous study j2]. Let us recall, for consistency, that 

^ we consider homogenization for the following elliptic problem 



- div ((^per(^) + &^(^, w)Cper(^))VtX,) = /(x) in P C M^ 

;e = on 9P, 

where the tensor Ap^r models a reference Z'^-periodic material which is randomly perturbed 
by the Z -periodic tensor Cper, the stochastic nature of the problem being encoded in 
the stationary ergodic scalar field bj^ (the latter getting small when r] vanishes). We have 
studied in [2| the case of a perturbation that has a Bernoulli law with parameter rj, meaning 
that brj is equal to 1 with probability 7] and with probability 1 — rj. In the present work, we 
address more general laws. The common setting is that all the perturbations we consider 
are, to some extent, rare events which, although rare, modify the homogenized properties 
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of the material. Our approach is a perturbative approach, and consists in approximating 
the stochastic homogenization problem for 

Aj^[X,UJj = Apgr[X) + bri [X,UJ) L^per 

using the periodic homogenization problem for Apg^- In short, let us say that our main 
contribution is to derive an expansion 

^; = a;,, + riA\ + r^^Al + o{v^), (1.2) 

where A* and A* are the homogenized tensors associated with A^^ and Aper respectively, 
and the first and second-order corrections A^ and A2 can be, loosely speaking, computed 
in terms of the microscopic properties of Ap^r and Cper and the statistics of second order 
of the random field 6^. The formulation is made precise in |2] and in Sectionsl2]andl3]below. 



Motivations behind setting (1.1), as well as a review of the mathematical literature on 
similar issues and a comprehensive bibliography, can be found in |2j. We complement our 
study of the perturbative approach introduced with [2] in two different directions. 

In Section |2] we rigorously establish an asymptotic expansion of the homogenized tensor 



in a mathematical setting where our input parameter (the field 6^ in (1.1)) enjoys appro- 
priate weak convergence properties, as r] vanishes, in a reflexive Banach space, namely 
a Lebesgue space L°^ (D , L^ (Q)) (with p > 1). In such a setting, we are in position to 
rigorously prove a first order asymptotic expansion (announced in |3] and precisely stated 
in [21 Theoreme 2.1] and Theorem [2] below) for the homogenization of Ajj, using simple 
functional analysis techniques very similar to those exposed in [1] . In our Corollaries p^ 
and|4] the expansion is pushed to second order under additional assumptions. 

Our aim in Section ^ is to further extend our formal theory of [2j . Recall that this 
formal theory, rather than manipulating the random field 6^ itself, consists in focusing on 
its law. We indeed assume that the image measure (the law) corresponding to the pertur- 
bation admits an expansion (see ( [3.4^ below) with respect to r/ in the sense of distributions. 
While [2] has only addressed the specific case of a Bernoulli law, we consider here more 
general laws and proceed with the same formal derivations. These derivations lead to a 



first-order correction ^^ in (1.2) obtained as the limit when A^ — )• 00 of a sequence of ten- 
sors A'l' computed on the supercell [— y, y] . It is the purpose of Proposition 7 to prove 
the convergence oi A^' . The second-order term A2 is likewise defined as a limit, up to ex- 
traction, of a sequence of tensors A2 when N ^ 00. The proof of the boundedness of the 
sequence A2' and thus of its convergence up to extraction is not given here for it involves 
long and technical computations. We refer to |T] for the details. As in |2], our approach in 
this Section exhibits close ties with classical defect-type theories used in solid state physics. 

We emphasize that, in sharp contrast to the exact stochastic homogenization of A„, the 



determination of the first and second-order terms in (1.2) relies on entirely deterministic 



computations, albeit of very different kind, for both approaches of Sections |2] and [3j 

Finally, a comprehensive series of numerical tests in Section |4] show, beyond those 
contained in p], that the two approaches exposed here are efficient and quite robust: 
the computational workload induced by the perturbative approach is light compared to 



the direct homogenization of [2j, and expansion (1.2) proves to be accurate for not so 



small perturbations. 

We complement the text by a long appendix. The reader less interested in theoretical 



issues can easily omit the reading of this appendix. Besides providing, in Sections 5.1 and 



5.2 and for consistency, some theoretical results useful in the body of the text, the purpose 



of this appendix is two-fold. We examine in details in Section 5.3 the one-dimensional 



setting, and we show that, expectedly, all our formal expansions can be made rigorous 



through explicit computations. We next demonstrate, in Section 5.4 that our two modes 
of derivation coincide in a particular setting appropriate for both the theoretical results 
of Section [2] and the formal results of Section [3j This final section therefore provides a 
proof of our formal manipulations of Section |3] in a setting - we concede it - that is not 
the setting the approach was designed to specifically address. Definite conclusions on the 
theoretical validity of the approach developped in Section [3] are yet to be obtained, even 
though applicability and efficiency are beyond doubt. 

Throughout this paper, and unless otherwise mentioned, C denotes a constant that 
depends at most on the ambient dimension d, and on the tensors Aper and Cper- We write 
(7(7) when C depends on 7 and possibly on d, Ap^r and Cper- The indices i and j denote 
indices in [1, dj. 

2 A model of a weakly randomly perturbed material 

For consistency, we first recall the general setting of our related work [2j. 

Throughout this article (il, T, P) denotes a probability space with P the probability 
measure and w G $7 an event. We denote by E(X) the expectation of a random variable 
X and Var{X) its variance. 



We assume that the group (Z , -|-) acts on and denote by r^, k E Z , the group action. 
We also assume that this action is measure-preserving, that is, 

V^ G J-,VA; G Z'^, P(^) = P(rfc^), 

and ergodic: 

V^ G F, (VA: G Z'^, ^ = TkA) =^ {F{A) = or P(^) = 1). 

We call F G L}^^{M.'^ , L^ {Cl)) stationary if 

\/k £ Ij ', F(x -\- k,uj) = F{x,TkUj) almost everywhere in x G M and u: £ Q. (2.1) 

Notice that if F is deterministic, the notion of stationarity used here reduces to 
Z'^-periodicity, that is, 

yk G Z'^, F{x + k) = F{x) almost everywhere in x G M''. (2.2) 

We then consider the tensor field from W^ x Q. to W^^'^: 

Ar^{x,Uj) = Aperix) + \{x,Uj)Cper{x), (2.3) 

where Ap^r and Cper are two deterministic Z'^-periodic tensor fields and 6^ a stationary 
ergodic scalar field. The matrix Aper models the reference periodic material, perturbed by 



Cper- This perturbation is random, thus the presence of 6,^. We refer the reader to [4] for 
a more detailed presentation of the stationary ergodic setting in a similar weakly random 
framework. 

We make the following assumptions on the random field hr^: 

3M>0,Vry>0,||6^||ioo(Q^f,)<M, (2.4) 

limiL°°(Q;L2(n)) -^ 0> (2-5) 

where Q is the unit cell [— g) 2]'^' 

Assumption (2.5) encodes that the perturbation for small ry is a rare event. Still, it is 
able to significantly modify the local structure of the material when it happens, for we do 
not require it to be small in L^(Q x fi) as ry — )• 0. 

We additionally assume that there exist < a < (3 such that for all ^ G M , for almost 
all X eW^ and for all s G [-M, M], 

a\C\^ < Aper{x)C ■ C, a\C\^ < {Aper + sCper) {X)C ■ C, (2.6) 



^per 



{X)^\ < m, I (Aper + sCper) {x)^\ < mi (2-7) 



We can therefore use the classical stochastic homogenization results (see for instance |7j 
for a comprehensive review or [2] for a concise presentation). The cell problems associated 



with (2.3) read 

' - div (AriiVw^ + a)) = in M'^, 
Vw^ stationary, E ( / Vw^ j = 0. 



(2.i 



Problem (2.8) has a solution unique up to the addition of a random constant. The 



function w^ is called the i-th corrector or cell solution. 



The homogenized tensor A* is given by 

yiell,4, A;e, = E(fAr,{Vw'^ + ei)Y (2.9) 

Throughout the rest of this paper we will denote by u;^ the i-th cell solution associated 
with Aper, defined up to an additive constant by 



(2-10) 



div {AperiVwf + a)) =0 in Q 
w^^ Z — periodic. 

The periodic homogenized tensor is then given by 



Vi G ll,dj, Al^,ei = [ AperiVw^ + e^). (2.11) 

JQ 

Due to the specific form of Afj, the following zero-order result can be easily proved. 
The proof is actually the same as that in Lemma 1 of f2], which relies on the fact that 
limiL°°(Q;L2(Q)) converges to as 7y tends to 0. 



Lemma 1. When r] —^ , A* —?' A. 



per- 



Our goal is to find an asymptotic expansion for A^ witli respect to rj, and a first answer 
is given by tlie following theorem announced as Tlieoreme 1 in [3] : 



Theorem 2 (Tlieoreme 1, [3j). Assume that hrj satisfies (2.4) and (2.5), and denote by 
"ir; = imilL°°(Q;L2(o))- There exists a subsequence of rj, still denoted rj for the sake of 
simplicity, such that ^ converges weakly-* in L^{Q; L'^{Q)) to a limit field denoted by bo 
when n — )• 0. Then 



for all i € [1, dj, the following expansion 



(2.12) 



holds weakly in L^{Q; L^(r2)), where w^ is the solution to the i-th periodic cell problem 
and v^ is solution to 



- div(^perV't;°) = div (6oCper(Vt(;° + Cj)) in R*^, 
Vv° stationary, E [ / Vf ° ) = 0. 

A* can be expanded up to first order as 

A*^ = Al^^ + mr,A\ + o{mr,), 
where 

Vi G \l,dl A\ei = [ E(bo)Cper{Vw^ + ei)+ f AperVE(7;0). 



(2.13) 



(2.14) 



(2.15) 



Proof. We fix i G [l, dj and define v. 



f] 



vj is solution to 



- div {A^Vv"/) = div (^Cper (Vw° + ei) ) in M^ 
Vvl stationary, E f / Vv^ ) = 0. 
Using an argument similar to that used in the proof of Lemma 1 in [2] , we have 

Vr/ > 0, W^vlWi^i^Qy^Q,) < -\\Cper (Vu;° + a) \\l^q)- 



(2.16) 



where a is defined in (2.6). 



The sequence Vf^ is bounded in L'^{Q x Q) and therefore, up to extraction, weakly 
converges in L?'{Q x O) to some limit which is necessarily a gradient and which we denote 
Vv^. Since 6^ converges strongly to in L^{Q x Q), brj^v^ converges to in T)'{Q x Q). 
It is then easy to pass to the limit r; — )• in (|2.16) and to deduce that vf is solution to 



' - div (AperVv^) = div {boCper {Vw^ + 6^)) in M"^, 

Vvi stationary, E f / Vv° ) = 0. 



Thus 



Vi»7-Vrt;V 
rrin 



converges, up to extraction, weakly to Vv^ in L'^{Q x Q). This amounts 



to say that we have the following first-order expansion: 

\7w''^ = Vt(;° + rrirjVvi + o{mr,) in L'^{Q x Q) weak. 



Inserting this expansion in (2.9), we obtain 



A*ei = A* ei + rrir, E{bo)Cper{Vw'i + ei)+mr, AperVEiv'D + oim^), 

JQ JQ 

which concludes the proof. 



D 



Remark 1. Notice that taking the expectation of both sides of (2.13), lE(wf) is actually the 
Ij'^ -periodic function that is the unique solution (up to an additive constant) to 



- div (^perVE(^O)) ^ ^i^ {E(bo)Cper {Vw^ + Ci)) in Q, 

E{v^) Z"^ - periodic. 



(2.17) 



The computation of A* up to the first order in m^ only requires solving 2d deterministic 



problems, namely (2.10) and (2.11), in the unit cell Q 



In fact, the situation is even more advantageous when Ap^r is a symmetric matrix, as 
shown by our next remark. 



Remark 2. Defining the adjoint problems to the cell problems (2.10), 



- div (^Jg^(Vu;° + ei)) = in Q, 
^Z;j Z — periodic. 



(2.18) 



where we have denoted by AlL^^ the transposed matrix of Aper, allows to write the first- 



order correction (2.15) in a slightly different form. Indeed, multiplying (2.11) by w^ and 
integrating by parts, we obtain 



[ AperVE{v^) • VZ^O = - / E(bo)Cper (V^° + 6^) • V 
JQ JQ 



^l 



Likewise, multiplying (2.18) by VE(u^) and integrating by parts yields 

[ AperVE{v^) ■ {Vw'j + ej) = 0. 
JQ 



Combining these equalities gives 



Aper^E{v^) ■ e, 



E(bo)Cper (Vw;° + Ci) ■ Vwl 



and thus (2.15) may be equivalently phrased as 

yii,j)ell,df, Alei-ej = f E{bo)Cper{Vw'! + ei)-iVw'j + ej). (2.19) 



When Aper is symmetric, w^ = Wp and solving the periodic cell problems (2.10) suffices 
to determine A* up to the first order in rrirj. 



Pushing expansion (2.14) to second order requires more information on b„: 



Corollary 3. Assume in addition to (2.4) and (2.5) that 

b^ = rjbo + r/2fo + o{ri'^) weakly -* in L°°(Q; L^{n)). 
Then 

• for all i G [1, d\, the following expansion 

"^wl = V-u;" + r/Vu° + rf'V zf + o{r]'^) 
holds weakly in L"^ {Q; L? [Q)) , where z? is solution to 

' - div(^perVz°) = div (foCper(Vii;° + 6^)) + div (5oCperVw°) in 
Vz° stationary, E f / Vzf ) = 0. 

• A* can be expanded up to second order as 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



where A\ is defined by (2.15) and for all i E [[l,(i]], 



^2^* 



Q 



nro)Cper{^w'i + ei) + r]' / Cp,,E(6oVt;y) + / ylperVE(zl'), (2.24) 



or equivalently, for all {i,j) £ [[l,(ip, 



■^2^i ' ^j 



E{ro)CperiVw^ + Ci) ■ {VWj^ + Cj) + / CperK(6oVwf ) • (VuJ/ + ej).(2.25) 



/ 



Proof. The proof follows the same patt ern a s that of Theorem [2] The computation of 
the second order relies on the fact that (2.20) implies that ^ converges strongly to bo in 
-L°°((5; L^(J7)), whereas the c onver gence was weak in Theorem^ Likewise, the expansion 

of the cell solution, namely (2.21), implies that — ^ ^ converges strongly to Vw^ in 

L2(g;L2(il)). We then obtain ( 12 23) ) and ( ]2.24[ ) by inserting ( |2.2l[ ) in (|2^, and deduce 

D 



( |2.25[ ) from \2.2A\ as in Remark 



i 



The computation of A* up to the order rf is much more intricate than that up to 
the order ry, for it requires determining E(feoVf ?). C omputing the periodic deterministic 
function E(z;^) solution to the simpler problem (2.17) is not sufficient in general. We have 



to determine the stationary random field v^ solution to (2.13) in 



It turns out that in a particular, practically relevant setting, we may still avoid solving 



the random problem (2.13). This setting presents the additional advantage to provide 



insight on the influence of spatial correlation. 

Corollary 4. Assume that brj is uniform in each cell ofL, and writes 

\{x,uj) = ^ lQ+fe(x)S^(rfcu;), 



(2.26) 



kei 



where B„ satisfies 



Vr?>0, llS^IUco(f^) <M, 
B^ = r]BQ + ifRo + o{7f) weakly in L2(ri). 



Assume also that 



Y^ \cov{Bo,Bo{Tk-))\ < 



oo. 



(2.27) 
(2.28) 



(2.29) 



Then the second- order term (2.25) can be rewritten 



A*e^ ■ Cj =E{Ro) / CperiVw"^ + Si) ■ (Vu5/ + Cj) + Var{Bo) / CperVti ■ {Vw] + Cj 



Q 



+ (EiBo))' / CperVs, ■ (Vw^j + Cj) 
JQ 

+ V COv{Bo, Boin-)) / CperVtii- - k) ■ {Vw^j + Bj), 



(2.30) 



'2 ^m>d\ 



where ti is a Lj^^{M ) function solving 

( - div {AperVti) = div {CpertQ (Vu;° + 6^)) in R'^ , 

[Vti£L\R'^), 



(2.31) 



and Si solves 



— div (AperVSi) = div [Cper (Vlfj + Cj)) in Q, 

Si TL — periodic. 



(2.32) 



defined in (2.20) here write 



Proof. We notice first that the specific form (2.26) of 6,j considered implies that bo and ro 

(2.33) 



bo{x,u)) = ^ lQ+kix)Bo{TkUj) 
kezd 

ro{x,uj) = ^ lQ+k{x)Ro{TkUj) 



(2.34) 



The rest of the proof mainly consists in showing that in this particular setting, Vv^ 
and the product boS/v^ can be written using the deterministic functions ti and Sj. The ex- 
istence of ti and its uniqueness up to an additive constant come from Lemmas 6 and 7 in [2] . 



We start by proving that the sum 



Y, {Boirkoj) - E{Bo)) Vti{x - k) 



(2.35) 



is a convergent series in L^{Q x 0). 



To this end, we compute the norm of the remainder of this series: 

2 



Y, {Bo{rk-)-nBo))Vtii--k) 



\k\>N 



L^iQxn) 



= E E cov{Bo{Tk-),Bo{Tr)) f Vt,(- - k)VU{- - I) 

\k\>N\l\>N ''^ 

\k\>N\l\>N 
^ Yl E l^°^(^0(■^'^■)'^0(^'•))lll^*^^- 



lL2(Q)) 



\lhq) 



\k\>N\l\>N 



< E l|Vt.(--A;)||i.(Q) ^ |co^(So(rfc-),^o(rr))| 

\k\>N' \ \1\>N 

< E l|V**(--^)lli2(Q) E \cov{Bo,BoiTi_k-))\ 

\k\>N \ \1\>N 

< E l|Vt.(--A;)||2,(Q) ^|co7;(5o,So(rfc.))|. 



Using (2.29), we obtain 



Y, {Bo{rk-)-nBo))Vti{--k) 

\k\>N 



<CY \\^tii--mhiQy (2.36) 



Since Vi, G L^(M'^), the right-hand side of (2.36) converges to zero when N goes to 
infinity. 



Consequently, (2.35) defines a vector T in L^{Q x Cl). It is clear from (2.35) that 
jf-^ = jpr" for all {n,p) G [1, dp. Thus T is a gradient, and there exists a function Vi such 



8Xn dXp 

that 



Vv^ = T + E(Bo)Vs, = Y (^o(rr) - lE(^o)) Vti(x - k) + E{Bo)Vsi. (2.37) 



Since Si is Z -periodic, we deduce from (2.37) that 



Vvi is stationary and E ( / Vvi I = 0. 



(2.38) 



We then compute, using ( |2.31| ) and (2.32), 

- diviAperVVi) = ^ -div(^perVti(--fc))(So(rfc-)-E(So)) 

kezd 

-diY{AperVSi)E{Bo) 

= Y div {CperlQ+k (Vu;° + ei)) (5o(rfc-) - E{Bo)) 



+dw{Cper{Vw^ + ei))E{Bo) 

Y div {CperlQ+k Bo{rk-) (Vu;° + e^)) . 



(2.39) 



Because of (2.33), ( |2.39 ) implies 

- d[v{AperVVi) = div {boCperiVw^ + Cj)) 



(2.40) 



It follows from (2.38) and (2.40) that Vi solves (2.13). As (2.13) has a solution unique 

(2.41) 



up to the addition of a random constant, we obtain 

Vvf = Vv^=Y^ (Boirk-) - lE(^o)) Vti(x - k) + E(^o)Vs,. 



We deduce from ( |2.33[ ) and ( |2.4lD that 
E(6oVt;°) = Y.Y1 ^Q+inBo{rrKBo{Tk-)-E{Bo)))Vti{- - k) 



+(E(5o))2 Y, 1q+;Vs, 



iezd 



J] Y, lQ+icov{Bo{n-),Bo{Tr))Vtii- - k) + {E{Bo)f ^ 1q+Ns,, 



and then that 



lQE(boVv^) =VariBo)Vti + J^ cot;(^o(-), ^o(rr))Vti(- - fc) 

fcezd,/c^o 
+ (E(^o))2Vsi. 



(2.42) 



We conclude by inserting ( |2.34[ ) and ( |2.42[ ) in ( |2.25[ ). 



D 

Theorem^ (and its two corollaries) are only of interest if E(6o) 7^ 0. Indeed, if E(6o) = 
it only states that A* = A* + o{mri). 

The prototypical case where Theorem [2] does not provide valuable information is the 
case studied in ||2]: 6^(x,a;) = XlfceZ'* ^Q+k{x)B^{u}) , where the B^ are independent iden- 
tically distributed variables that have Bernoulli law with parameter ?], i.e are equal to 1 
with probability r] and to with probability 1 — r/. Then, using the notation of Theorem [2| 
6^ = 6,,, nirf = y^ and 6o = 0, and we only get A* = A*g^ + o(y^) (while appendix 6.1 
of p] shows that there exists a tensor A\ such that A* = A* + r]A\ + o{ri) at least in 
dimension one). Omitting the dependence on the space variables since 6^ is uniform in 
each cell of Z in this particular setting, a suitable functional space F on fi to obtain a non 
trivial weak limit of inpfj— would be L^{Q) for the norm of each B^ in L^{^) is equal to 
r]. The Dunford-Petti weak compactness criterion in that space is however not satisfied by 

converges in the set of bounded measures 



The reason is of course that 



ii"»)iiii(n) IINIlLi(n) 

to a Dirac mass. The techniques used in the proof of Theorem |2] and its two corollaries 

thus do not work in this setting. 

The above considerations somehow suggest that an alternative viewpoint might be 



useful. Because of (2.5), the image measure dP!S of ^(x,-) converges to a Dirac mass in 



the sense of distributions. Our alternate approach, related to our work [2j, consists in 
working out an expansion of the image measure (or of the law) , rather than an expansion 
of the random variable. Like in |2], our manipulations are mostly formal. Some rigorous 
foundations, in specific settings, are provided in the appendix. 



10 



3 A formal approach 



3.1 A new assumption on the image measure 

For simplicity, we assume as in Corollary UJ that brj is uniform in each cell of Z , and is of 
the form 



b^{x,uj)= Y, T^Q+kix)B'^{uj), 



(3.1) 



k& 



where the B^ are independent identically distributed random variables, the distribution of 
which is given by a "mother variable" Brj. For convenience we slightly modify (2.27) and 
require 



3e>0,Vr/>0,||5^||ioo(o)<M-e 



IB 



»?llL2(n) -^ . 



(3.2) 
(3.3) 



Assumption (3.2) is a technical assumption which implies in particular that for every 



T] > 0, the image measure dP^j of -B^ is a distribution with compact support contained in the 
open set ] —M, M[. Of course the specific values of M and e have no particular significance. 
Throughout the sequel we denote by £'(] — M,M[) the space of distributions on M with 
compact support in ] — M, M[, and by (T, c^) the action of a distribution T £ £'(]— M, M[) 
on a test function tp G C°°{] — M, M[). Basic elements of distribution theory are recalled in 



Section 5.1 of the appendix, for convenience of the reader not familiar with technical issues. 



that for every ^ £ C°°{] - M,M[), 



Because of assumption (|3.3|) and Lebesgue dominated convergence theorem, it is clear 



Since K{ip{Brj)) = {dPrj,ip) and (/^(O) = {So,ip) where 5o is the Dirac mass at 0, dPr^ con- 
verges to 6o in £'{] - M,M[). 

This leads us to assume that dPr^ satisfies 

dPr^ = 6o + r]dPi+r]'^dP2 + o{if) m£'{]-M,M[), (3.4) 

which is equivalent to 

yip G C°°(] - M, M[), E{ip{B^)) = {dP^, if) = ip{0) + vidPi,^) + V^{dP2, v) + o{7]^). 

Of course dPi and dP2 also have a compact support contained in ] — M, M[ : for every 
test function ip with compact support in M\[— M + e, M — e], it holds for all ?] > 

{dPrt, ^) = n^{Br,)) = = ri{dPi,^) + rf{dP2, ^) + o{rf), 

which yields (dPi, y?) = {dP2, f) = 0. Then the supports of dPi and dP2 are contained in 

[-M + e,M-e] c]-M,M[. 

Denoting by M' = M — e/2, we deduce from Proposition Il3^ of the appendix that 
there exists a constant C > and integers pi and p2 (namely the orders of dPi and dP2 
respectively) such that 



yipeC°^{]-M,M[), \{dPi,^)\<C sup sup 

se[-M',M']0<n<pi 



ds^ 



fis] 



(3.5) 
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Vv3 G C"^(]-M,M[), |(dP2,^)| <C sup sup 

xe[-M',M']0<n<p2 



ds-^ 



Lp{s) 



(3.6) 



Let us now give some additional motivations underlying assumption (3.4). 



The first motivation is related to our work presented in [2] in which B^^ has Bernoulli 
law with parameter r/, meaning that it is equal to 1 with probability r] and with prob- 
ability 1 — r/. Then the image measure dP^ is equal to 5q + r]{5i — (5o), so that it satisfies 
(3.4) exactly at order 1 with dPi = 5i — 5q. 



The second motivation comes from the following result, which shows that there is an 



easy way, used in our numerical experiments, to build perturbations satisfying (3.4). 



Lemma 5. Consider B a random variable in L^iyt). Let K he a positive real, and define 



B.q = 7/i?l|,^5|<^. Then B^j, which obviously satisfies (3.2) and (3.3), also satisfies (3.4) 
with 



dPrj = 6o- 7]E{B)6'o + '^E{B^)6'^ + 0(7]^) in £'{R). 



(3.7) 



Proof. Let us denote by dP the image measure of B, and consider ip G P(M) (i.e ip G C°^(M) 
and has compact support). Then 



(3.8) 
(3.9) 



{dP^,(p) = / ip{r]s)dP + ^{0) / dP 



<f{'ns)dP+ I {^{0) - (f{r)s))dP. 

\r)s\>K 



Since B is in L^(0), 



/ dP = 0{v'), 

J\ris\>K 

and thus, ip being a bounded function, 

{dPr„^)= l\{r^s)dP + 0{rf). 
Then, since if £ 2?(]R), there exists C > such that 



VsG 



Again using B G L^(ri), this implies that 



^{^s)-^{0)-WiO)-\s'p"iO) 



<C7]\sf. 



V(r?s) - (^(0) - 7?V(0) - ^sV(O)^ 



dP ^0 



which is just a rewriting of (3.7) since / dP = 1, / sdP = K{B) and / s dP = K{B ). D 

Before exposing our approach in this new setting, we prove the following elementary 
result which we will often use in the sequel: 
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Lemma 6. It holds {dPi, 1) = and {dP2, 1) = 0. 

Proof. It holds on the one hand (dP^, 1) = 1 since dP^i is a probabihty measure, and on 
the other hand 



{dPr,, 1) 



so that the conclusion follows. 



{6o, 1) + vidPi, 1) + vHdP2, 1) + o(r?') 
l + ri{dPi,l)+r,\dP2,l) + o{r,^), 



D 



3.2 An ergodic approximation of the homogenized tensor 

Let us consider a specific realization u) G Q of Ajj in In = [~"2") y] i ^ being for simplicity 
an odd integer, and solve the following "supercell" problem: 

- div (Arj{x,u}){'\7w^' ''^ + ei)]=0 inliy, 
wl^'^ {NZf - periodic. 
Then we have 



(3.10) 



\fi ell, dj, A;e, = Jim^^Efj Ar,{x,oj){Vw'^''''''{x) + e,)\ dx. (3.11) 



The proof of (3.11) is given in [2]. We only outline it here for convenience. We know 



from Theorem 1 in [5j that 



1 



N'i Ji 



ri,N,u] 



Arf{x, cj)(Vu;^' '^{x) + Cj) dx converges to A*ej almost surely in w G Jl. (3.12) 



Since 



Nd 



, v,^,^ 



In 



Afj{x,oj){'Vw^' '"(x) + ei)dx is the periodic homogenization of A^(x,a;) on 



In, it is also well known that for all {i,j) G [[l,(i] , 

j^i A-\x,oj)dxj e^-Cj < j^ Ar,{x,u){Vw'^''^''^{x) + a) ■ Cjdx 



< 



^ ( / Ar^{x,oJ)dx] a -ej. 



so that for all N G 2N + 1, for all r/ > and for almost all a) G il. 



1 



r),N,u], 



—^ Aj^{x,uj){Vwf ' {x) + ei)-ejdx 

'J In 



</3, 



(3.14) 



we can take the expectation in (3.12) and get (3.11). 



where /3 is defined by (2.7). Using (3.14) and the Lebesgue dominated convergence theorem. 



Remark 3. The same result holds for homogeneous Dirichlet and Neumann boundary con- 
ditions instead of periodic conditions in the definition of w^' '^ (see ^5] for more details). 

For convenience, we label the unit cells of In from 1 to A^ . The k-th cell is denoted 
by Qk, for 1 < k < N . A given realization Ari{x,ub) can then be rewritten 



Ar^{x,Uj) = Aper{x) + ^ Ig Jx)SfcCper(a:), 



fc=l 
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with Sfc = B^{uj) for all k G [1, A^ J. The B^[oj) being independent random variables, the 



Af" 



joint probability of the A^ -uplet (si, • • • , Sjqd) is simply the product I I dPri{sk) 



k=l 



Remark 4. The approach exposed in the sequel works also, with minor changes, for random 
variables which are not independent but correlated with a finite length of correlation. We 
present it in the independent setting for simplicity. 



Afd 



We now define A'^^- ^'n^ = Aper + Y^ Iq^SkCper for (si, • • • , s^d) G [-M, M]^\ We 

fe=i 
denote hy w^ ' ' ^ the solution of the z-th cell problem for the periodic homogenization 
of A^^'"'''^N<i on Ijy, that is 



l.-,Sjvd 



divL4*i'-'"iv'^(Vw;/' '^"+ei) =0 Iu/tv, 



w. 



si,-- ,s^d 



{NZf - periodic. 



Then, defining 



we have 



■^r} ^i 



Nd 



E 



An{x,io){Vw'^' '"^(x) + Cj) I dx, 






(3.15) 



(3.16) 



^L(//^^'-'"^(^-^^^^'""^^^^On/^^(-)- (^■^^) 



It is proved in Lemma 14 of the Appendix that 'Vw^ ' ' ^ is a C°° function of (si , • • • , s^d) 
in ] — M, M[ . Thus, since dPi and dP2 have compact support in ] — M, M[ (as well as <5o 



of course), we can make these distributions act on A^'^'"'''^Nd and ^w- 

of (Sl,--- ,Sj^d). 



as functions 



It follows from (3.4) that 

jyd j^d j\fd jyd 

lldPr,isk)=ll6oisk)+vY.dPiisi) n ^oisk) 

fc=l k=l 1=1 k=l,k=/=l 

+ Y E E dPl{si)dPi{Sm) n ^0{Sk) 

1=1 m=l k=l,ky^{Lm.} 

^d ^d 

+ n^J2dP2{si) n So{sk)+ON{v^) m£'{]-M,M[^''). 

1=1 k=l,k^l 



(3.18) 



We stress that the remainder oi\f{r] ) in (3.18) depends on A^, hence the notation. 



Moreover the products (3.18) are to be understood as tensorized products: we work in 

£'{] - M, M[) 01 £'{] - M, M[) 02 • • • ®m-i £'{] - M, M[) C £'{] - M, Mf'). 



Inserting (3.18) in (3.17), we obtain the following second-order expansion 



\*,N 



,N 



\*,N , 2 A*,N 



Al'"+r^Ar+v'Ar+ON{v'). 
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(3.19) 



Before making the first three orders in (3.19) precise, note that (3.11), (3.16) and (3.19) 
imply 



,*,N 



,N , 2 A*,N 



A*= hm U^'^+r/^r +r/Mr +ojv(r/^ 



(3.20) 



In the sequel we exchange in (3.20) the limit in A^ and the series in rj in order to guess 
a second-order expansion of A* depending only on ry. Since we are not able to justify this 
permutation, our approach is formal. 



We now detail the first three orders in (3.19) 
First, we notice that for i G [l, rf]]. 



A*^^e- 



\fc=l ■'^t^ 

^ In 



w, " + ei) 



J In 
= 4* p- 

which obviously gives the zero-order term expected for A* Then 



A^ Cj 



7V° / N°- \ 

i_i \ 7 1 i.-Li J In I 



1=1 \ k=l,k=/=l 

It is easy to see that, by (A^Z) -periodicity of w^ 

'In 



i.-.s„ 



dPiisi) n ^o{sk), A^^'■■■''^N''iVw■"^^^'''''+e^) 

I 1 :.-Li J In I 



k=l,k=/=l 



does not depend on I. The expression (3.21) can then be rewritten 



,*,N 






We change the notations for convenience, and define, for s G [—M,M], 



j^sfl _ ^s,0---,0 



-Aper ~T~ 'S JJ- Q L^per ) 



and w- = w- solution to 



(3.22) 



(3.23) 



,5,0/^7 i-,s,0,N 



divL4t''^(V<*'"'"+e.) =0 in In, 



w. 



l,s,0,N 



{NZy - periodic. 



(3.24) 



The matrix A^' corresponds to the periodic material with a defect of amplitude s lo- 



l,s,0,N 



cated in Q (i.e at a position G Z in Ijv), and w^' ' ' is the i-th cell solution for the 
periodic homogenization of A^' in /jv- Since w-'^' ' = w^' ' ' , it is of course a C°° 
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function of s g] - M,M[. 

With these notations, we find that 



,N 






AY'e, = {dPi{s),J^ AfiVwl'^'^^'^' + e, 



(3.25) 



For the second-order term, we first define the set 



iV-1 iV-1 



r^ = |fcGZ^Q + A;C/7v} 



The cardinal of T/v is of course A^"^, and I J {Q + k} = In- 



(3.26) 



feeTjv 



For {s,t) G [-M,M]'^ and k eTN.'^e define 



A. 



s,t,0,k 



Aper ~r SJioCper + t^Q+k^peri 



(3.27) 



and w- solution to 



w. 



iS,t,0,k 



tfl,k^yj^^2,s,tfl,k,N 

I 

2,s,t,0,k,N ^j^T^\d 



divUf'"''^(Vt/;f™"+e.. 



in/TV; 



(3.28) 



(NZ)" - periodic. 



The matrix A2 ' ' corresponds to the periodic material with two defects of amplitude 
s and t located in Q and Q + k (i.e at positions G Z and fc G Z in I^) respectively. 
The function w^ '*'''■ is the z-th cell solution for the periodic homogenization of A2 ' ' 
in In- It is a C°° function of (s,t) g] - M,M[^. 



Then computations similar to that presented for the first order yield 



A.2 S-i 



\ Y. (dPi(.)ciA(t), / Af°'^Vu;^^'*'°''='^ + e.) 



k£TN,k¥=0 



+ (dP2{s),j^ Af(V^MAA^ + e,) 



(3.29) 



A setting with zero, one and two defects is shown in Figure [T] in the two-dimensional 
case of a reference material Ap^r consisting of a periodic lattice of circular inclusions. 







•••# •••• 

















Figure 1: From left to right: zero defect, one defect and two defects. 
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Remark 5. It is illustrative to consider the particular case where the random variable B„ 



has a Bernoulli law. This is the case treated in jSj/. Then, expansion (3.4) holds exactly 
with dPi = 6i — So- The distribution dP2 and all other terms of higher order identically 



vanish. The expressions (3.25) and (3.29) then coincide with (3.17) and (3.18) in /^. 



,N 



In the next section we prove that A-^ converges to a finite hmit when A^ — )• oo. The 



.A' 



case of the second-order term A2 , which is shown to be a bounded sequence and thus to 



converge up to extraction, is discussed in Section 3.4 



3.3 Convergence of the first-order term 



,N 



We study here the convergence as N goes to infinity of A^ defined by (3.25). 



.*,N 



Proposition 7. The sequence A-^ converges in M to a finite limit A\ when N — )• 00. 



*,N 



Proof. We fix {i,j) G [[l,d]] and study the convergence of A*' a ■ Cj. 



Using (3.24) and tlie adjoint problems defined by (2.18), we first obtain, for all 



se [-M,M], 



+ ei) ■ e,- 






,0,N 



+ e^)- (ei + Vw°). 



Then, letting the distribution dPi act on the left and right-hand sides, and using (3.25), 
we find that 

Al'^'ci ■ e, = /dPi(s), I AfiVwl''''^'"' + ei) ■ (e,- + Vu>°)\ . (3.30) 

Because of the definition of A^' , 

/ Af{Vwl''''''' + e^)■{ej + Vw^)= f ApeAVw]'''''''' + e,) ■ {e, +Vw 

J In J In 



+ / sCper(Vt«l'^'°'^ + ei) ■ (e,- + V^2;°). 
JQ 



(3.31) 



Next, using (2.18), 



/ A,e.(V^.l'^'°'^ + e.) • (e, + Vu;°) = / (Vu;,^'^'"'^ + e.) • A^Ae, + Vw^ 
Jlisr JIn 

= I ei-Af;,,iej + Vw^). 
J In 

We know from Lemma p^ that {dPi, 1) = 0. Thus 



(3.32) 



dPi{s), I ei-yl^,^(ej+Vu;°) 



Collecting ( |3.30[ ), ( |3.3lD , ( |3.32D and ( |3.33[ ), we get 



At{" ei ■ Cj = (dPi{s), j sCper{Vwj 



+ ei) • {cj + Vw° 



We now define 



l,s,0,N l,s,0,N 

Qi =Wi"' -Wi. 



(3.33) 



(3.34) 



(3.35) 
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l,s,0,N 1 
q^ solves 



div [Al'^Vql'"'^' ) = div(slQCper(Vw° + ej)) in In, 



^i,sfl,N ^j^^^d _ periodic. 



(3.36) 



Using (3.35) in (3.34), we rewrite 



A^ Ej • ej 



sdPi(s), / Cper(Vli;° + Ci) • (Cj + Vtt>°) 






l,s,0,Af 



,7,0^ 



+ ( dPi(s), \ sCperiVql'''''''' + a) ■ (e, + V*," 



(3.37) 



The rest of the proof consists in showing that 






(dP,{s), J sCperi^q^'''''^'' + a) ■ (e,- + Vwl 



which is of course equal to 



l,s,0,Af 



sdPi{s), I Cper(Vg/'^'^'^^ + ei) ■ {ej + Vw'^) 

converges to a finite limit when A^ — )■ oo. 

More precisely, defining 

Vs e [-M, M], ViV G 2N + 1, f^{s)= I CpeAVql'''"^'^ + e^) • (e^- + Vw'j), 

JQ 

we will prove that the sequence f^ and its derivatives converge uniformly, when A^ goes 
to infinity, to a limit function /°° and its derivatives. 



Applying Lemma 15 of the appendix to (3.36), we obtain that for all s G [— M, Ml, 

l,sfi,N ^ - ' - 

■ converi 

function solving 



VQj''^' ' converges in L'^{Q), when A^ — )• oo, to Vq-'^' '°°, where g^''^' '°° is a Lf^^(R.'^) 



- div (^f Vg^^'°'°" 



divislgCperiVw^ + a)) in M^ 



(3.38) 



Moreover, arguing as in the proof of Lemma [15] (given in our previous work |2j), it is 
easy to see that for all A; € N and all s G [—M,M], VSg'g^''^' ' converges in L'^{Q) to 

VT Ofc 1,5,0,00 

We then define /°° by 

Vs G [-M, M], r{s) = f Cpe.(V(?^^'°'°° + e,) • (e, + V*"). 

Jo 



Because of (5.4) and (5.5) in Lemma 16 of the appendix, and using a classical result of 



differentiation under the integral sign, it is clear that 






[ C7per.(V5f=(?^^'°'^ + e,) • (e, + Vw^), 
Jo 
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and 



VA: 



e nys G] - M, M[, ^r{s) = f C7per.(V5,^(?^^'°'°° + e,) • (e, + Vwf). 

J CJ 



l,s,0,oo 



The convergence of VS^g^''^' ' to Vdgq-'^' '°° in L'^{Q) for every A; G N thus yields 



VA; G N,V5 G] - M,M[, Mm -^f{s) = -r^r{s). 

N-^+oo as'' as" 



(3.39) 



On the other hand, we deduce from Lemma 17 that there exists a constant C{pi,M) 
(recall that pi is the order of dPi{s)) such that for all k G [[0,pil|, 

V(s,/)G]-M,M[2,ViVG2N + l, |— /^(,) - — /^(s')| < C(pi,M)|. - s'|. (3.40) 



It is straightforward to see that (3.39) and (3.40) imply that 



d^ ^ 

VO < k < pi, — rf converges uniformly to -—rf°° in 1 — M, M[. 
ds'^ 



d'' 
ds^' 



(3.41) 



It follows from (3.5) and (3.41) that 

{sdPi{s),f{s))^{sdPr{s),r{s)), 
and then 









l,s,0,N 



( dPi(s), / sCper(Vg,^'^''^'" + ei) ■ (e, + Vw') 



N^oc 



l,s,0,oo 



dPiis), / sCperiVqr'°" + ei) ■ (e, + V^^^ 



(3.42) 



Collecting (3.37) and (3.42), we conclude that A[' converges to a limit tensor A[ 
defined by 



V(i,i) G [1, df, Ala ■ e, = /sdPiis), j Gper (Vu;° + e^) • (e, + Vw^j] 



{' 



l,s,0,oo 



+ ( dP,{s), I sCper ( Vg,^'^'^'°° + e, ) • (e, + Vu;°) ) . 



(3.43) 
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3.4 Second-order term 



,A^ 



For completeness, we state here the result concerning the second-order term A2 in (3.19 ), 
proved in jT]: 



Proposition 8. The sequence A*2 defined by (3.29) is bounded in W^^'^ and therefore 
converges up to extraction. 



:,Ar 



We firmly believe that A2 is actually a convergent sequence, as shown by our numer- 
ical tests thereafter. We also stress that the explicit computations of Section [5. 3| prove this 
convergence in dimension one. 
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4 Numerical experiments 

The purpose of this section is to assess the numerical relevance of the approaches of Sections 
[2]and|3] To this end we build and homogenize stochastic composite materials using laws 
that satisfy the assumptions of these sections. Our motivations are not strictly identical for 
the two approaches. In contrast to the first approach which relies on a rigorous proof, our 
second approach is formal and we thus need to demonstrate its correctness experimentally 
(note that the tests performed in [2] in the Bernoulli case are already to be considered as 
a component of the validation of the approach). We wish to check that the expansions 
derived in Sections [2] and |3] provide an accurate and efficient approximation to the direct 
stochastic computation. The limited computational facilities we have access to impose 
that we restrict ourselves to the two-dimensional setting. We first explain our general 
methodology, which is the same as that presented in |2|, and then make precise the specific 
settings. 

4.1 Methodology 

We mainly consider as in [2] a reference material Ap^r that consists of a constant background 
reinforced by a periodic lattice of circular inclusions, that is 

Aperixi,X2) = 20 X Id + 100 ^ lB(fc^o.3)(2;i,a;2) x Id, 

where 5(^,0.3) is the ball of center k and radius 1. Loosely speaking, the role of the 
perturbation is to randomly eliminate some fibers: 

Cper{xi,X2) = -100 ^ ^BI,k,Q.Z){xi,X2) X Id. 



We will also, in our last test, consider a laminate 

Aper{xi,X2) = 5 + 10 ^ lKa;i<i+l(xi, X2), 

with the perturbation yielding an error in the lamination direction: 

Cper{xi,X2) = 10 ^ \l<x2<l+l{xi,X2) X Id - 10 ^ l/<a;,</+i(xi, X2) X Id. 

For both materials (shown in Figure ^ , we have chosen the values of the coefficients in 
order to have a high contrast between Aper and Ap^r + Cper and thus for the perturbation 
to have an important impact on the microscopic structure. The specific value of these 
coefficients has no other significance. 



We will consider different perturbations 5.^, all of which satisfy (3.1 ) with the B^ inde- 
pendent and identically distributed. 

Our goal is to compare A* with its approximation A* + rjA'l' + ff-^ ■ A major 
computational difficulty is the computation of the "exact" matrix A* given by formula 



(2.9). It ideally requires to solve the stochastic cell problems (2.8) on M . To this end we 



first use ergodicity and formula (3.11 ), and actually compute, for a given realization lo and 
a domain /^r chosen here to be [0, N]'^ for convenience, A^' ioj) defined by 

K;''i^y^ = ^J A,{x,uj){Vwl''^-{x) + e,)dx. (4.1) 
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Figure 2: Left: a periodic lattice of circular inclusions. Right: a one-dimensional laminate. 

In a second step, we take averages over the realizations oj. 

For each uj, we use the finite element software FreeFem++ (available at www.freefem.org) 



to solve the boundary value problems (3.10) and compute the integrals (4.1). We work 



with standard PI finite elements on a triangular mesh such that there are 10 degrees of 
freedom on each edge of the unit cell Q. 

We define an approximate value Ar^' as the average of A^' (w) over 40 realizations uj. 
Our numerical experiments indeed show that the number 40 is sufficiently large for the 
convergence of the Monte-Carlo computation. We then let A'^ grow from 5 to 80 by steps 
of 5. We observe that A^' stabilizes at a fixed value around A^ = 80 and thus take j4^' 
as the reference value for A* in our subsequent tests. 

The next step is to compute the zero-order term A*^^ and the ffist-order and second- 
order deterministic corrections. Using the same mesh and finite elements as for our refer- 
ence computation above, we compute A* using (2.10) and (2.11). The computation of 
the next orders depends on the setting: 



in the setting of Section 2, the ffist-order correction is given by (2.15) in Theorem [2] 
and is thus independent of A'^; since h^q is of the form (2.26), we use formula (2.30) 



in Corollary |4] for the second-order correction which depends on N through the term 
ti defined on M by (2.31), and which has to be approximated on Ijv; we let A^ grow 



from 5 to 80 by steps of 5; 



*,N 



,*,N 



in the setting of Section 3, the corrections A^^' and A2 are respectively given by 



(3.25 ) and ( 3.29[ ); we let A^ grow from 5 to 80 by steps of 5 for j4^' ; the computation 



of A2 being far more expensive (there is not only an integral over /^v but also a sum 



over the A^ cells in (3.29)), we have to limit ourselves to A'^ = 25 and approximate 



the value for A'^ larger than 25 by the value obtained for A'^ = 25. 
We stress that there are three distinct sources of error in these computations: 
• the finite elements discretization error; 
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the truncation error due to the replacement of M with In, in the computation of 



the stochastic cell problems (2.8) that are replaced with (3.10), as well as in the 
computation of the integrals (4.1 ); 

• the stochastic error arising from the approximation of the expectation value by an 
empirical mean. 

Detailed comments on these various errors and the way we deal with them are provided 
in [2j. We just emphasize, in the setting of Section^ that it is not our purpose to prove 
through our tests that 

a; = a;,, + rjAl + rfA% + o{if) 

with a 0(77^) which would be independent of A^, of the number of realizations and of the 
size of the mesh. We only wish to demonstrate that the second-order expansion is an ap- 
proximation to A* sufficiently good for all practical purposes. We will observe that A*2 
is not only bounded as stated in Proposition ^ but converges to a limit ^2, and that both 
yl*' and A*2 converge to their respective limits faster than A!^' to A* (which is expected 
since the former quantities are deterministic and contain less information). We will also 
observe that A* -|- rjA*^ is closer to A* than A* and that the inclusion of the second 

order improves the situation for A*^ + ^^1' + ^^^2' ^^ even closer. 

To present our numerical results, we choose the first diagonal entry (1,1) of all the 
matrices considered. Other coefficients in the matrices behave qualitatively similarly. We 
illustrate a practical interval of confidence for our Monte-Carlo computation of A* by 

showing, for each iV, the minimum and maximum values of ^^' (w) achieved over the 40 
realizations u. 

We will use the following legend in the graphs: 

• periodic: gives the value of the periodic homogenized tensor A*^.] 

• first- order: gives the value of the first-order expansion; 

• second-order: gives the value of the second-order expansion; 

• stochastic mean, minima and maxima: respectively give the values of A!^' and the 
extrema obtained in the computation of the empirical mean. 



Finally, the results are given for various values of ry which serve the purpose of testing 
our approach in a diversity of situations, and in particular for perturbations that are "not 
so small". 

4.2 An example of setting for our theory in Section [2] (and |3]) 

Consider B.^ = 7/G1o<7jG<i where G is a normalized centered Gaussian random variable. 
It is easy to check that 

Br^ = 7]G1q<g<+oo + o{rj^) in L^{Vt), 
so that Corollary [4] of Section 2 applies. Alternatively, we can use Lemma [5] which gives 

dP, = <5o-r?-^<5^ + ^5^' + o(7?2)in£:'(M), 
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to perform our formal approach. We verify in Section 5.4 of the appendix that both ap- 
proaches yield the same results up to second order. 

We show results for the lattice of inclusions and for r] = 0.1 and rj = 0.2 (Figures p] and 
E] respectively) . 

The results are very satisfying for both values of rj. The first-order correction, which 
does not depend on N, enables to get substantially closer to A* Moreover, it is clear 

(especially from the close-ups) that the second-order correction A2 converges very fast 
(convergence is already reached at A^ = 5), and in particular much faster than the stochastic 
computation A^' . It also provides excellent accuracy. 

4.3 A first example of setting for our formal approach of Section [3] 

Consider i?^ a random variable having Bernoulli law with parameter rj, and G a normalized 
centered Gaussian random variable independent of Rrj. We define the product random 
variable B^i = Rr) x rjGluiQ\<i. Then 

E{cp{B^)) = E(99(i?^ X 7?Gl|^,G|<i)) 

2 
= r?(^(0) + r/E(G)(/.'(0) + ^^"(0) + o{v')) + (1 - r?)^(0) 

= ^(0) + |^c^"(0)+o(r?3). 



This implies 



3 

dP^ = 5o + ^5o + o(r/3) inf'(]R). (4.2) 



In this case we only consider the first-order correction since the dominant order in (4.2 ) 
is already tiny. We present the results in the case of the lattice of inclusions, for r/ = 0.2, 
r] = 0.3 and r] = 0.5 (Figures n] m u^ respectively) . 



Once again, our approach converges rapidly and allows for an accurate approximate 
value of A* even for r] as large as 0.5. 
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Figure 3: Inclusions - results for a Gaussian perturbation and t] = 0.1. Above: complete 
results. Below: close-up on Arj' and the first and second-order corrections. 
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Figure 4: Inclusions - Results for a Gaussian perturbation and rj = 0.2. Above: complete 
results. Below: close-up on Arj' and the first and second-order corrections. 
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Figure 5: Inclusions - results for perturbation (4.2) and r] = 0.1. Above: complete results. 
Below: close-up on A^' and the first-order correction. 
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Figure 6: Inclusions - results for perturbation (4.2) and r] = 0.3. Above: complete results. 
Below: close-up on A^' and the first-order correction. 
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Inclusions - Bernoulli times Gaussian - eta=0.5 
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Figure 7: Inclusions - results for perturbation (4.2) and r] = 0.5. Above: complete results. 
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Below: close-up on Ajj' and the first-order correction. 
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4.4 A second example of setting for our formal approach of Section [3] 

Consider Rrj a random variable having Bernoulli law with parameter rj, and U a uniform 
variable on [0, 1] independent of i?^. We define B^ = R-q — f]U . Then 

n^{Br,)) = E{ip{Rr, - rjU)) 

= r]E{ip{l - r]U)) + {1 - i])E{^{-rjU)) 
= r? (99(1) - r?E(C/)v9'(l) + 0(7?)) 

+(1 - rj) ((^(0) - r?E([/)(^'(0) + ^E(C/2)^"(0) + oi^. 

= ^{o) + v{-HU)ip'{o) + ip{i)-ip{o)) 

so that 

+ 7?2 (-E{U){5[ - 5',) + \nU^)5"{Q)^ + oiv') in ^'(M). ^^'^^ 

Notice that this complex case is a mixture of Sections 2 and 3. The first-order per- 
turbation is of course only the sum of the first-order perturbations for a Bernoulli law 
(Section 3 and |2]) and a uniform law (Section 2). The interaction of these laws at order 2, 
and notably the S[ term, is much more involved and requires the computation of the cross 
derivatives of w- '^' ' ' ' with respect to s and f at s = and t = 1. 

We give the results in the case of the inclusions and for rj = 0.05, rj = 0.1 and rj = 0.2 



(Figures pi ^ 10 respectively). 



For rj = 0.05 and rj = 0.1, the results display the same features as in our previous tests 
and are very good. The case rj = 0.2 is instructive: the second-order expansion signifi- 
cantly departs from the "exact" value provided by the direct stochastic computation. Our 
interpretation is that, far from contradicting the validity of our expansion in the limit of 
small rj, it shows the limitations of the approach. The value r] = 0.2 is too large for the 
expansion to be accurate in the case of a lattice of inclusions with a high contrast between 
the inclusions and the surrounding phase. 

Interestingly, a value of rj twice as large (0.4) provides a very accurate approximation 



for another material, as shown by our final test performed on the laminate (Figure 11). 



Our approach has limitations and deteriorates, like any asymptotic approach, for large 
values of rj. The threshold is case dependent. The approach is however generically robust. 
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Figure 8: Inclusions- results for perturbation (4.3) and rj = 0.05. Above: complete results. 
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Below: close-up on Arj' and the first and second-order corrections. 
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Figure 9: Inclusions - results for perturbation (4.3) and rj = 0.1. Above: complete results. 
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Below: close-up on ^^' and the first and second-order corrections. 
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Figure 10: Inclusions - results for perturbation (4.3) and rj = 0.2. Above: complete results. 
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Below: close-up on Arj' and the first and second-order corrections. 
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Figure 11: Laminate - results for perturbation (4.3) and t] = 0.4. Above: complete results. 
Below: close-up on A^' and the first and second-order corrections. 
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5 Appendix 

The objectives of this appendix are diverse. We first quickly recall some elements of 
distribution theory. We then prove technical results used in Section [3j Next we show that 
the approach formally derived in Section [3] is rigorous in dimension one. Finally we prove 
that this approach is also rigorous, in general dimensions, in a specific setting close to that 
of Theorem [2] and Corollary |4j 

5.1 Elements of distribution theory 

We recall here some basic definitions and results of distribution theory for convenience of 
the reader. See |6] for a comprehensive presentation. 

In this section O denotes an open set in M. 

Definition 9. We denote by D{0) the space of infinitely differ entiahle functions on O 
having compact support in O. 

Definition 10. T is a distribution on O if T is a linear form on T>{0) satisfying the 
following continuity property: for every compact K C O, there exists an integer p and a 
constant C such that for all Lp G 'D(0) having compact support in K, 



\{T,^)\<C sup 

x^Kfi<n<p 






(5.1) 



The space of distributions on O is denoted by T>'{0). 



If the integer p in (5.1) can be chosen independently of K, the distribution T is said to 



have a finite order. The smallest possible value for p is called the order of T . 

Definition 11. A distribution T G T)'{0) is said to have compact support if there exists a 
compact set K C O such that for all ip £ T>{0) having compact support in 0\K, (T, (/?) = 0. 

The support of T is defined as the smallest compact set K which satisfies the above 
assertion. 

The space of distributions on O having compact support is denoted by £'{0). 

Proposition 12. If T £ £'{0), its action on V{0) can be naturally extended to C°°{0). 
Denoting by K a compact neighborhood of the support ofT, and by x o- cut-off function in 
T>{0) equal to 1 on the support of T and vanishing on 0\K, we define 

Vv9gc°°(0), (r,vp):=(r,x(/^). 

This definition does not depend on K and x- 

Proposition 13. // a distribution T is in £'{0), it has a finite order. Denoting by p its 
order and by K a compact neighborhood of the support of T, there exists a constant C > 
such that: 

'ipeC'^iO), |(T,v?)| <C sup '^" 

x&K,0<n<p 



dx" 
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5.2 Some technical results 

This section is devoted to the proof of technical lemmas used in Section[3] Loosely speaking, 
these lemmas all deal with the variation of the supercell correctors defined by (3.15), (3.24), 



and (3.28) with respect to the amplitudes of the defects. 



Lemma 14. Let Hp^j.{Ij^) be the set of {N'Z)'^ -periodic functions in Hj^^^iJR.'^) with zero 
mean on Ij^. The function 



where w- 






£ H (Im), 



w- 



Ju 



w- 



'^ and w- 



is defined by (3.15), is C°° . 



Proof For (si, • • • , s^d) G [-M, M]^ , w^ 



-Sir--, 



is the unique solution to 



div(^^i'-'"ivd(Vu;f''"''^'+ei; 



w"''"'''^^{NZ)'^- periodic, 



w. 



in /jv , 



0, 



In 



SO that F is well defined. 

Let us now define G :] - M, ufxHl^^ilN) -^ H-^{In) by 



G(si,--- ,sj^d,w) 

si,---,Sj^d 



-div {A^^'"' '^N-i (\/w + a)) , 
so that F{si, • • • , s^d) = w^ ' ' ^ is the unique solution to 

G(si, • • • , Sj^d,F{si, • • • , s^d)) = 0. 
It is easy to see that G is a C^ function, and that 

V/i G ^p^,,(/iv), d^G{su--- ,Sj^d,w)-h = -diyiA'''-''N'iVh), 
where dwG{si, ■ ■ ■ , Sj^d,w) is the first derivative of G with respect to ti; at (si, • • • , Sj^d,w). 

The Lax-Milgram theorem and the coercivity of A*i'"''^jv'' show that dwG{si, • • • , S]sfd,w) 
is an isomorphism. We can therefore apply the inverse function theorem and deduce that 
-F is C^, with dsiF the unique solution to 

- div (A'^'- ^'N'i^Vds.F)) = div {iQ^CperiVF + e^)) iniN, 

ds^F (NZ)'^ - periodic, I ds^F = 0. 

Jin 

Arguing by induction, we obtain that F is a C°° function. 



n 



For consistency, we state next a lemma proved in [2l Lemma 6] 



Lemma 15. Consider f G Lp'{Q), and a tensor field A from W^ to W^^'^ such that there 
exist A > and A > such that 

Ve G M'^, a.e in x G M'^, A|^|2 < A{x)i ■ C and |^(x)^| < A|^|. 

Consider q solution to 

r - div {AVq"") = div(lQ/) in In, 
\ q^ (NZy - periodic. 
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Then Ij^Vq converges in L (M ), when N goes to infinity, to Vg°°, where q°° is a Lj^^(Ji 
function solving 



- div {AVq'^) = div(lQ/) in R'^, 



(5.3) 



Lemma 16. Consider q^'^' ' and q-'^' '°° solutions to (3.36) and (3.38) respectively, and 
A; G N. There exists a constant C{k,M), such that 

VsG]-M,M[, ViVG2N + l, WVO^ql'^'^^'^'h^^j^) < C{k,M)\\Vw^ + ei\\L2^Q), (5.4) 



Vs e]-M,M[, \\Vd^ql'''''^h2^^,) < C(A.,M)||Vw;° + e.|U2(Q). 



(5.5) 



Proof. Multiplying the first line of (3.36) by q^ ' ' ' and integrating by parts, we find that 



l,sfl,N 



llVg. 



l,s,0,Af| 



\lHIn) 



< M- 



\a 



per\\L°°{Q) 

a 



Vu;° + ei||i2(Q), 



(5.6) 



where a is defined by (2.6). 



Thus (5.4) is true for A; = with C(0, M) = M 



\c, 



per||L°o(Q) 



1,S,0,N 



Next, the first derivative dsq^' ' ' is solution to 

- diviAfVdsql'"'^''') = div {igCperiVw^ + Ci)) + div (igCperVql'''^'^ 

dsql'''^'^ {NZ)'^ - periodic, 
from which we deduce 



in I, 



N' 



(5.7) 



\^dsq, 



l,s,0,N< 



and, using (5.6), 



IL^Im) 



l,s,0,N, 



< 



\a 



per\\L°°{Q) 

a 



K-'n 



l,s,0,Af| 



V7 U I II I IIVT i,S,U,JV II 

Vwi +ei||i2(Q) + llVg^' 11^2 



(Q) 



IV7Q J-,s,u,iv|| ^ W^PerWL^^iQ) ,,. iMiv7 , || 

|V9,g.' WmiM) ^ (^ + l)l|Vu;i +ei||i2(Q). 



a 



Thus (5.4) is true for fc = 1 with C(l, M) = {M + 1 



\\Cper\\ L°° (Q) 



Finally, we have for k >2 

- diviAfVd'^ql'''''''') = kdiv ( iQCperVar ^g, 



,fc-l^l,s,0,Ar 
per ^ ^s 



in / 



N, 



Vd^ql'''^'^ {NZ)'^ - periodic, 



so that an easy induction proves (5.4). The proof of (5.5) is identical. 



(5.8) 



n 



The following result is an immediate consequence of Lemma 16 



Lemma 17. Consider q^'^' ' and q^''^' '°° solutions to (3.36) and (3.38) respectively. For 
every /c G N, there exists a constant C{k, M) such that for all (s, s') g] — M, Mp, 

ViV G 2N + 1, ||V5,^g^^'°'^ - V5,^g^^''°'^||i2(,^) < Cik, M)\\Vw^ + e,h2^Q^\s - .'|,(5.9) 
WVdUl''''^^ - Vd'J'''''n\LHu-) < C{k, M)\\Vwf + e.b2(Q)|. - .'|. (5.10) 
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5.3 The one-dimensional case 



We address here the one-dimensional context. All the computations are explicit, for the 
settings of Sections 2 and [3] To stress the fact that we deal with scalar quantities, we 
use lower-case letters for the tensors. Note also that in this section Q = [—\,V\ and 

5.3.1 An extension of Theorem [2] 

The following theorem extends the result of Theorem^ stated in L°°{Q] L^(ri)), 
to L°°(Q; LP{yi)) for any p g]1, oo]: 



Theorem 18 (one-dimensional setting). Assume that d = 1, that 6^ satisfies (2.4) and 
v ■~ limiLoo(r-i ^-LPifi)) ~^ ^ f^''^ some p > 1. There exists a subsequence of rj, still 



m. 



rj^O 



denoted r/ for simplicity, such that ^ converges weakly-* in L°°([— g, 5]) -^^(^)) to a limit 
field denoted by 69 when r/ — t- 0. Then 

• the expansion 



d 



d 



-w ' 



„o 



d 



,0 



w + m„——v + o{m 



VJ 



(5.11) 



dx dx ' ^ dx 

holds weakly in L^([— ^, ^];LP{Q)), where w^ is the periodic corrector and v'^ solves 



d d Q 



d f 7 / d 



d n I f^ d r, 

-—V stationary, E / -—v 
dx \J~- ^^ 



0. 



m . 



(5.12) 



a* reads 



^ „o , 1^ , _ P . d^,.. 



* = a* + m^ / "&{bo)cper{-^w +l) + mr^ I aper-r-^{v ) + o{mr,). 



'dx 



dx 



Proof The periodic and stochastic correctors can be computed explicitly. They are re- 
spectively given by 



A. 

dx 



w 






""per 



I - 1 and 



A. 

dx 



w' 



E 



-1 



Note that w^ is in H^^'°°(-i, \). 



We define v'^ 



-. It solves 



m^ 



d , d „. d fb„ , d r. . 



dx 



v^ stationary, E 



1 dx 

2 



0. 



m 



(5.13) 
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We deduce from (5.13) that 



d 



mr. 



^per\ 



d 

dx 



w^ + l) + kn, 



(5.14) 



where A:^ depends only on w. Since /c^ is by construction stationary ergodic, it is constant, 
and we compute from (5.13) and (5.14): 



1 



E 



Wlr. 



1 

2 1 



X IE ''" ^^-^ ^".,,0 



i ar, 



1 a 



L-per 



^ 



<i^'"" + ')l- 



Since w^ is in W^^°^' ^ ^ 



2 ' / \ " 2 

2, 2)1 fl»7 is coercive and Cp^r is bounded, it holds 

lifer. 



^■r)\ 



< c- 

< c. 



^'7iiLi([-ii]xn) 



m, 



7? 



l^'?llLi([-i|]xfi) 



5r?llL-([-i,|];LP(n)) 



This implies that A;„ is a bounded function of 77 whatever p > 1 and thus, using (5.14), 



that ^v^ is bounded in L'^{[-\, \]]LP{Q.)) for all p > I. As a result, for p > 1 
converges weakly and up to extraction in L^([— j, j]? -^^(^)) to a limit we denote -^v^ 



~^,v 






1 ii 



The random field 6^ tends to in L ([— ^, ^]; ^^(rj)). Since it is bounded in L°°([- 2, 2J '^ 
Q), it converges to in L?'{[—i^, 2]i-^'^(^)) foi' ^-^ r > p. By Holder inequality it also 
converges to in L^([— 2, 2]) -^'^(^)) foi' all 1 < r < p. Thus it converges to in in 



1 1 

'2' 2 



|;L5([7))where'gi^. 



The space L^Q.i^ i]. ^g(j7)) being the dual of L2([-i, i]; LP(1))), we obtain that 
brjCper ^V^ tends to in V{[—^,l]x Q,). We can then take the limit ry — )• in (5.13) 
and obtain that v^ is solution to (5.12). 

We have thus proved that ^^ i.dx'^^ ~ dx^^^ converges, up to extraction, weakly to 
^v^ in L^([-5, \];U'{Q)), which is equivalent to (|5.1l|). 



dx 



The second assertion of Theorem 18 is obtained by inserting (5.11 ) into the expression 

^^ of a;. n 

Note that the proof of Theorem [18] depends crucially on the fact that we are able to 
solve explicitly the cell problems. 

Theorem 18 allows for a better intuitive understanding of Theorem [2] In dimension 
one, the homogenized coefficient is explicitly given by 



E 



— 1 Oiper ~r Or; Cper 



which, when br^{x,Lo) = > lr^^^^]^i(x)-B^(rfca;), may be rewritten as the formal series 



fcez 



f2i-i)'mB,)') 



»? k=0 



^per 



jiper 



per' 



(5.15) 
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Assume now that there exists p > 1 such that ||-B„|| rp(n\ — J- when ?7 — t- and n „ n '' 

f II //lll^'-|,S2j / \\-Dr]\\LP{n) 

converges weakly in U'{n) to some Bq with ¥j{Bq) ^ 0. We have in particular 



\B. 



r?llLP(Q) 



E(So) / 0, 



which, since E(|i3^|P) — t- 0, implies 

n\Br,?)= o (E(B,)). (5.16) 

r;-s>0+ 

We now claim that, without loss of generality and up to an extraction in r/, we may take 



P 

E(|i?^|2l 



2 in (5.16). Indeed, if p < 2, then since Brj is bounded in L°°(r2), (5.16) implies 



sequence |rg 



o (E(i3„)). On the other hand, if p > 2, we consider the normalized 

r?-s>0+ 

in Lp'ip.). Up to extraction, it weakly converges to B2 G L^(0). Since 



Brt 



n\\L'2{n) 



nsr,) 



nBr,) 



\B. 



»7llL2(f7) 



IB 



»?llLp(n)) 



|Bry||L2(n) l|-Sr;||LP(n) 



where the left hand side converges to E(i3o) 7^ and 



\B. 



\\Br,\\ 



'jll£2(Q) 



E/' R 

inequality, E(52) = lim t— 

v^o \\B, 



We then take p = 2. Since E(|B^p) 



LP(n) 



is bounded by 1 by Holder's 



7^ and (5.16) is satisfied with p = 2. 



n\B,\'') 



o (E{Br,)) and B^ is bounded in L°°(0), 



o (E(5^)) for all k>2. 



This intuitively expresses that all orders higher than or equal to 2 are negligible as 



compared to the first-order term in the series (5.15), and thus that a kind of "separation 



of scales" is satisfied. This is of course formal since one has to check that the remainder 
term consisting of the sum of all terms of order higher than or equal to 2 is o (E(i?^)), so 
that 



+ 



1 a. 

1 



per 



1 a 



per 



-E(Bj-^)+o(E(i?J) 



a, 



per 



-1 



-2 



+ m. 



1 CLper 



1 ar. 



E(5o 



^per 



''per 



+ 0{¥.{Br,)). 



But this is the purpose of the proofs of Theorems [2] and 18, using another viewpoint, to 
show this is indeed the case. 



5.3.2 The setting of Section [3] in dimension one 

We now prove that our approach of Section [3] is rigorous in dimension one. 
Lemma 19. In dimension d = 1, it holds 

a* = a* + 7]a[ + rfa^ + o{rf), 



where a\ and 03 are the limits as N ^ 00 of a^' and a^' defined generally by (3.25) and 
(3.29) respectively. 
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Proof. Recall that in dimension one, a* is given by the simple explicit expression 



E 



-1 



A ^T)GT I (J/fiLji 



dPr,{s), 



per 1^ "ri'-per 



-1 



— Cirpgy I (bC-^ 



■per 



The proof thus consists in inserting expansion (3.4) in this explicit expression and iden- 



tifying successively the first three dominant orders. 



K) 



Using (3.4|), we write 



1 



1 a, 



+ r,idPi{s), 



■per 



— 1 Ciper ~r SCp^f 



+ lf{dP2{s), 



+o{rf) 



= {0*per) M ^ + ^"*per ( dPl{s), 

+o{rf). 
This yields the expansion 



A ^T)C.T i oO' 



per 



+ rj^a*{ dP2{s) 



1 Qjper ~r SCp^f 



A ^T)GT i oC"] 



-per 



a^ '^per VK^per) \ ^^ly^J'i 



1 (XxfCT \ SCij 



+ r]\a;,,f{dP,{s), 



r]'{a*r(dP2{s), 



"pe 



1 Q'per ~r SCp^j- 

1 

5 1 



1 Qjper ~r SCp^j' 



(5.17) 



+ o(r/^ 



We now devote the rest of the proof to verifying that the coefficients of 77 and r/^ in 
(5.17) are indeed obtained as the limit as A'' — )• 00 of a*^ and Cg' defined generally by 



(3.25) and (3.29) respectively, in this particular one-dimensional setting. 

The function w^''^'^'^ generally defined by (3.24) satisfies here 
d 



sflf d l,s,0,N 
dx V°^ ^^^^ 

w, iV — periodic. 



+ 1) 



in 



N N 



(5.18) 



We easily compute using (5.18): 
.s,0/ d 



a 



^ dx 



_yjl,Sfi,N 



+ 1) 



N 



N ^ 

T 1 

JV Gper + sliri uCper 



N{N{a;,,)-'-f{s)) ' 



"per 



/-I I 

"'per ' 



''perf J./ \ , V^per) 



f{sf + o{N'^), 



where f{s) 



sc 



per 



i dperyCiper i SCper) 
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Thus a{ defined generally by (3.25) takes here the form 



*,N 



2 / 



= Na;,, {dP.is), 1) + (a;,,)2 (dPi(s), lis)) + o(l). 
We know from Lemma 6 that (^dPi{s), l) = 0, whence 

"*''' ' al = {a;,,f{dPi{s),fis)). (5.19) 



N~^oo 



N-1 Af-li 



Likewise, we compute from (3.28), for A; G [- 2 ; 2 Ji' 

/ — 
-'<tfl,kf d 2,s,tfi,k,N -- ^^ I f '^ 1 



ar-(^^w'"" 



+ 1) 



N 



-1 



N aper + sl[_l_l]Cper + il[fc_ l ^fc+ij'^per ^ 



* \-l 



sc 



■per 



N{N{a;,,)~^-f{s)-f{t))~\ 



tc, 



■per 



(XperXP'per ~r tCpef) 



Then 



s,t,0,k, « ^,2,s.t,0.k,N 
"2 I . ^^ 



* -^3 



Notice that this expression is i ndep endent of k (and so of the distance between the two 



,N 



defects), so that Oj' defined by (3.29) here reads 

,M N{N-1) 



dP,{s)dPr(t),a;,,. + ^^^{f{s) + fit)) + ^^(J{s) + fit))' 



2 V ---v^/--^v/'--per . ^ 

+ N(dP2is),a*+^-^^fis))+oil). 



(5.20) 



to 



Since we know from Lemma 6 that (^dPiis), l) = and (^dP2is), l) = 0, (5.20) reduces 



*,N _ ( * \3 



Qr, = a 



per / 



dP,is),fis)f + (a;,,)2 <rfp2(,)^ j(,)> + 0(1). 



Thus 



,A^ 



Af-s>cx) 



(a;,j3 (^dP,is),fis)f + (a;,,)2 {dP2is),fis)) . 



(5.21) 



Finally, since /(s) 
we have 



2 1 



1 a 



per 



1 Oper ~r SCpgr 



{dPiis)Jis)) = -{dPiis), 



and (dPi(s),l) = {dP2is),l) = 0, 



(5.22) 



and 



{dP2is)Jis)) 



dP2is) 



1 tlper ~r SCp^f 



1 dper \ SCper 



(5.23) 



In view of ( [5Tf| , ( [539| , (|5^, (|5^ and ([5^, we have proved 

a* = a*g^ + 770^ + r]'^a*2 + 0(77^). 



n 
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5.4 A proof of the approach of Section [3] in a specific setting 

The purpose of this final section is to prove that the formal approach of Section |3] is rigor- 
ous in a setting related to that of Corollary |4j 

More precisely, we assume that the random field br^ satisfies the assumptions of Corol- 
lary [4] These assumptions do not imply that the image measure dPrj satisfies assumption 
(]3.4[) which is at the heart of the approach of Section ^ so that we have to impose that 



dPr) additionally satisfies (3.4). The following preliminary result then gives the necessary 



form of the expansion of the image measure dP^. 
Lemma 20. Assume that 6,j satisfies 

br,{x,uj) = Y, tQ+kix)B^iu;), (5.24) 



where the B^ are i.i.d random variables, the distribution of which is given by a "mother 
variable" B^ satisfying 

V7?>0,||S^|Uoo(f^)<M, (5.25) 

Br, = 'i]Bo + r]^Ro + oiv'^) weakly in L^ ((7). (5.26) 



Assume further that the image measure dPr, of B„ satisfies (3. 4)- Then 



2 
dPrj = 60- vnBo)6'o + yE(B2)<5^' - r]^E{Ro)6'o + 0(7?') in £\R). (5.27) 



Proof. Firstly, notice that -^ converges strongly to Bq in L (fi) because of (5.26). Now 
consider ip € P(M). We have on the one hand 

and on the other hand 

E {B^^iBr,)) = r]{s^dPuip) + if{s^dP2, if) + o(r?2). 

Thus s'^dPi = and s^dP2 =E{B^) 5o in P'(M). It is then well known that there exist 
71, Ki, 72, K2 ill IR such that 



E(m 



0; x'f 
0- 



dPi = 71(^0 + '^I'^o ^"^^ dP2 = 72(^0 + «:2'^o "I :^ — "^1 

Lemma |6] implies 71 = 72 = 0. Then, we have 

E(BJ = r]E{Bo) + r/2E(i?o) + o(??') 

and also 

E{Br,) = visdPi, 1) + v^{sdP2, 1) + 0(772). 

Thus {sdPi, 1) = E{Bo) and {sdP2, 1) = E (flo), from which we deduce ki = -E(^o) 
and K2 = — E(^o)- 

D 
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Theorem [2] and Corollary |4] rigorously yield the second-order expansion 

A*^ = a;,, + r]Al + 7]^ A* + o(7?2) 



with Al and A^ respectively defined by (2.15) and (2.24). 

On the other hand, using (5.27), Section l3| yields the formal expansion 



A*^ = a;,, + rjAl + r,'A* + o{r,') 



,*,N 



where A^ is the limit of the sequence A^' defined by (3.25) or equivalently by (3.30), and 



*,N 



A^ the limit (proved only up to extraction) of the sequence A2 defined by (3.29) 



The rest of this section is devoted to verifying that A\ coincides with A\ and A2 



coincides with A2 in the specific setting of Lemma 
5.4.1 First-order term 



20 



Using (5.27), (3.30) reads 



A*fe, • e, = -E(5o) / <5^(s), j 5Cper(Vu;^^'°'^ + e,) • (e, + Vu;^) 
and we compute 

A*'^ = E(^o) / Cper(Vii;^°'°'^ + a) ■ (e, + Vw^A. 
JQ 

Setting s = in (3.24), it is clear that w-' ' ' is equal to the periodic corrector wf. Then 

(5.28) 



,N 



Al'" = E{Bo) / CperiVw'i + a) ■ (e, + Vu;^). 



,N 



Clearly A-^ does not depend on A^ and its limit is then 



Al = E(5o) / Cper(V^y + e,) • (e, + Vw),^). 



(5.29) 



We recognize in the right-hand side of (5.29) the first-order coefficient in (2.19), which 



we know from Remark ^ is equivalent to (2.15). Theorem k] therefore shows that the 
first-order expansion 

a; = Al^, + riAl + 0(77) 

is correct with the values of the coefficients given by our formal approach of Section [3] 

We now proceed similarly with the second-order coefficient. 
5.4.2 Second-order term 



Using the adjoint cell problems (2.18) in (3.29) as in the proof of Proposition uj let us first 
rewrite 



-^2 ^* ' ^i 



Y, (dP,{s)dPi{t), f sCperVt.;,''^'*'"'"'^ • (e, + Vu;") 

T., u^n \ -J Q 



keTN,k^O 



l,s,0,N 



+ ( dP2{s), / sCperiVw^'"'''''' + Ci) ■ {c, + Vw^) 



(5.30) 
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Inserting (5.27) in (5.30), we start by focusing on 






( dP2{s), / sCperiVwl'"'"''' + 6,) • {cj + Vri)^ 



l,s,0,N 






l,0,0,Af 



l,s,0,N 



Denoting by dgW- ' ' ' , the first derivative of w,' ' ' evaluated at s = 0, we compute 



l,s,0,Af 



dP2is), / sCperiVw^'''''''' + €,) • {Cj + Vw 



E(S2) / CperVdsw]'''''''^ ■ {Vw'j + e,) +EiRo) I CperiVwl'^'^'^ + a) ■ {Vw'^ + ej] 



1,0,0,N /yy ~0 



E{B^) / CperVdsW^'"^''''' ■ (Vu;^ + ej) + E{Ro) / Cper{Vw^ + e.i) ■ {Vw^^ + e^). 



It follows from (3.24) that dsW,' ' ' solves 



i,o,o,Af 



- div{AperVdsWi' ' ' ) = div [iqCperiyw^ + Bj)) iniN, 

dswl'^'^'^ {NZf - periodic. 



(5.31) 



A^ -;> oo, to Vti defined by (2.31 



Applying Lemma |15| to (5.31), we deduce that VdsW^^' ' ' converges in L'^{Q), when 

in Corollary PI Consequently, 



Q 



l,sfl,N 



dP2{s), / sCper I Vw- ' ' ' + Cj ) • (Cj + Vw 



„r,0^ 



^ / Af->oo 



E(52) / CperVti ■ (Vu;° + ej) + E(iJo) / Cper {Vw^ + ei) ■ (V*° + Cj). (5.32) 



Next, we address 



, 2,s,t,0,k,N , , V7 ~0 



Denoting by Sjw; 
at s = t = 0, we have 



^ /dPiis)dPi{t), I sCperVl 

E(5o)2 ^ ((^o(^)<5o(i)> / sCperViof ^•*'°'''^ • (e, + Vu;0) 

with respect to t evaluated 



the first derivative of w^- 



Y, (dPi{s)dPi{t), I sCperVw^'-'''^'"''''' ■ {e, + Vwl 
mo? V / C7,e.V9,u;f °'°'°'^'^ • (e, + V^"). 



keTN,k^O 



(5.33) 



k£TN,k^O 



It follows from (3.28) that 9tt«' ' ' ' ' solves 



div{AperVdtW^' ' ' ' ' ) = div {lQ+kCper{VWi + Si)) in/AT, 



dtw. 



2,0fl,0,k,N 



{NZ? - periodic. 



(5.34) 
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Defining df = Y^k^TN dtWi' ' " ' , 
that solves 



it is easy to see tliat d^ is a Z -periodic function 



- diY{AperVdf) = div {CperiVwf + d)) in Q, 

dfZ^- periodic. 



(5.35) 



wliere Si is defined by (2.32) in Corollary 4 



Since problem (5.35) has a unique solution up to an additive constant, Vd] = Vs 



l,0,0,Ar 



Finally, comparing (5.31) to (5.34) for A; = 0, we find that Vdtw 



2,0fl,0fl,N 



is equal to 



VdgWj^' ' ' and then also converges in L {Q) to Vij when N — ?■ oo. 



Then, starting from (5.33), 



J^ /dPi{s)dPi{t), I sCperVu;^ ''*'°'^'^ • (ej + Vw 



keTN,k^O 



(E(5o))' [ Cper V Va, 



2,0,0,0,fc,Af 



keTN 



(E(So))' / Cpe.V9t«;, 
JQ 



2,0,0,0,0,7V 



(e,+V*0) 



(cj + Vw] 



: (E(So))' / CperVsi • (e, + V^°) - (E(5o))' / CperV5iu;f °'°'°'°''^ • (e,- + Vu;°) 

^ (E(^o))^ / CperVSi • (ej + VU)°) - (E(So))^ / CperVti ■ (cj + V^°) (5.36) 

-^°° JQ JQ 



,N 



It entails from (5.30), (5.32) and (5.36) that A2 converges to a limit A2 defined by 

A*ei ■ ej = E(i?o) / Cper (Vu;° + e^) • (V-w^ + ej) + Var{Bo) [ CperVU ■ (Vu;° + ej) 
JQ ^Q 

+ (E(^o))' f Vsi-ie.+Vw'l). 
JQ 



A2 is equal to the second-order term given by (2.30) in Corollary since we deal with 
independent random variables in each cell of Z . Thus the second-order expansion 

a; = a;,, + r,Ai + ifA*2 + o(7?2) 

derived from the formal approach of Section |3] is correct in this specific setting. 
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